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Abstract 

The crystals of potassium hydrogen carbonate (KHCO3) and the KDCO3 analogue are isomorphous. 
They are composed of hydrogen or deuterium bonded centrosymmetric dimers (HCO3 )2 or (DCO^^- 
The space group symmetry of KH p Di_ p C03 (p w 0.75) determined with neutron diffraction is identical 
to those of KHCO3 and KDCO3. This is at variance with a random distribution of H and D nuclei. These 
crystals are macroscopic quantum systems in which protons or/and deuterons merge into macroscopic 
states. 
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Introduction 



The formalism of quantum mechanics extrapolated from the level of elec- 
trons and atoms to that of everyday life leads to conclusions totally alien 
to commonsense, such as Schrodinger's Cat in a superposition of "alive- 
dead" states and nonlocal observables [TJ [2J. Such conflicts lead to a 
dichotomy of interpretation consisting in that, while at the microscopic 
level a quantum superposition indicates a lack of definiteness of outcome, 
at the macroscopic level a similar superposition can be interpreted as sim- 
ply a measure of the probability of one outcome or the other, one of which 
is definitely realized for each measurement of the ensemble [21 IH [5] ■ For 
open systems, this dichotomy can be legitimated by decoherence [5] stip- 
ulating that an initial superposition state should lose its ability to exhibit 
quantum interferences via interaction with the environment. However, 
since the quantum theory does not predict any definite dividing line be- 
tween quantal and classical regimes, macroscopic quantum behaviour is 
possible for systems decoupled from, or very weakly coupled to, the sur- 
roundings [7J. In principle, there is no upper limit in size, complexity, and 
temperature, beyond which such systems should be doomed to classicality. 

As a matter of fact, defect- free crystals are macroscopic quantum sys- 
tems with discrete phonon states at any temperature below melting or 
decomposition. This is an unavoidable consequence of the translational 
invariance of the lattice. Yet, the dichotomy of interpretation arises for 
O-H- • • O hydrogen bonds, when the coexistence of two configurations, 
say 01 — H • • • 02 and 01 • ■ • H — 02, is conceived of as "disorder" [8]. 
Protons are thought of as dimensionless particles, with definite positions 
and momenta, moving in a double-well coupled to an incoherent thermal 
bath. By contrast, vibrational spectra provide unquestionable evidences 
of phonon states showing that the translational invariance and the quan- 
tum nature of lattice dynamics are not destroyed by proton transfer. It 
is therefore compulsory to elaborate a purely quantum rationale avoiding 
any classical ingredient [9|. 

Such a theoretical framework has been elaborated for the crystal of 
potassium hydrogen carbonate (KHCO3) made of hydrogen bonded cen- 
trosymmetric dimers (HCOj - ^ separated by K + entities [TOlfTT] . The co- 
existence of two configurations for protons was first regarded as statistical 
in nature [121 113] , However, vibrational spectra (infrared, Raman, inelas- 
tic neutron scattering-INS) and neutron diffraction evidence the macro- 
scopic quantum behaviour of protons and theoretical arguments suggest 
that this behaviour is intrinsic to the crystal state [TH Q~51 [TBI [IZl HE] • 

In the present paper, we examine whether macroscopic states survive 
in a mixed crystal KH p Di_ p C03, grown from a water solution composed 
of H2O, HDO, D2O molecules. Since there is no segregation of H and D 
atoms into separate domains or clusters, it might seem straightforward to 
suppose that the crystal is a mixture of HH, HD, DD dimers, with proba- 
bilities p 2 , 2p(l — p), (1 — p) 2 , respectively [T9l [20] . A random distribution 
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of such dimers should destroy the translational invariance of the crystal, 
as well as — needless to say — the space group symmetry, and this should 
be evidenced by neutron diffraction. 

We report below neutron diffraction measurements showing that a 
mixed crystal (p w 0.75) is isomorphous to the pure analogues KHCO3 
and KDCO3. This is at variance with any random distribution of H and 
D nuclei, so the concept of macroscopic states previously elaborated for 
pure crystals should hold for isotope mixtures, as well. 

The organization of this paper is as follows. The crystal structure of 
KHCO3 at 300 K is presented in Sec. Q]and the theoretical framework for 
H or D macroscopic states is presented in Sec. [5] In Sec. [3] we present 
evidences of macroscopic proton states obtained with neutron diffraction. 
Finally, some consequences of the P2\/a space group symmetry deter- 
mined with neutron diffraction for the mixed crystal are discussed in Sec. 

H 

1 The crystal structure of KHCO3 and KDCO3 




Figure I: Schematic view of the crystalline structure of KHCO3 at 300 K. 
The arrows point to the sites occupied at low temperature. Dashed lines 
through protons are guides for the eyes. The ellipsoids represent 50% of 
the probability density for nuclei. 

The crystal of KHCO3 is monoclinic, space group P2i/a (CiyJ, with 
four equivalent entities per unit cell (Fig. [T]). Centrosymmetric dimers 
(HCC>3~)2 linked by moderately strong hydrogen bonds, with lengths Roo ~ 
2.58 A, are well separated by K + ions. All dimers lie practically in (103) 
planes, hydrogen bonds are virtually parallel to each other, and all pro- 
tons are crystallographically equivalent (indistinguishable). This crystal 
structure is unique to probing proton dynamics along directions x, y, z, 
parallel to the stretching (fOH), the in- plane bending (<50H), and the 
out-of-plane bending (7OH) coordinates, respectively. 

From 14 K to 300 K, there is no structural phase transition. The 
increase of the unit cell dimensions and of the hydrogen bond length are 
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marginal. Only the population of proton sites changes significantly. Below 
rs 150 K, all dirners are in a unique configuration, say L (see arrows in 
Fig. [T]). At elevated temperatures, protons are progressively transferred 
along the hydrogen bonds to less favored sites (configuration R) at « 0.6 
A from the main position. The center of symmetry is preserved and all 
proton sites remain indistinguishable. The population of the less favored 
site (or interconversion degree p ~ 0.22) is determined by an asymmetric 
double- well [3TJ [53] , but there are controversies as to whether this leads 
to statistical disorder [T5J [T3] or quantum derealization [TTJ [T51 [231 [5] • 

The isomorphous P2\/a crystal of KDCO3 is made of a very similar 
arrangement of centrosymmetric dimers. The hydrogen bond length is 
slightly longer (i?oo ~ 2.61 A) and the population of the less favored 
site is smaller (p w 0.12) [IT] . This suggests a slightly greater potential 
asymmetry for KDCO3 than for KHCO3 [21] . 

2 Theory 

In this section we present the bases of the theoretical framework leading 
to macroscopic states for protons or deuterons [T7] [IS] . 

2.1 The adiabatic separation 

Within the framework of the Born-Oppcnhcimer approximation, the vi- 
brational Hamiltonian can be partitioned as 

= Wh/D + Wat + C( H /D)at, (1) 

where 7Yh/d an d Hat represent the sublattices of light nuclei (H + or D + ) 
and heavy atoms, respectively, while C(H/D)at couples the subsystems. For 
OHO or ODO bonds, coupling terms between OH/OD and O- • • O degrees 
of freedom are rather large [5] [H] , hence beyond the framework of the 
perturbation theory. Two approaches, either semiclassical or quantum, 
are envisaged. 

Semiclassical protons or deuterons are dimcnsionlcss particles moving 
across a potential hypersurface [25l [26l [27] [28] [29] . Complex trajectories 
involving heavy atom coordinates lead to mass renormalization, and to 
incoherent phonon-assisted tunnelling [T2J, [T3l [30] . This approach is quite 
natural when the Born-Oppenheimer surface is known, but quantum ef- 
fects can be severely underestimated. 

Alternatively, if the classical concept of "trajectory", totally alien to 
quantum mechanics, is abandoned, adiabatic separation of the two sub- 
systems, namely Hh/d an d W a t, may lead to tractable models [51 [Tfl |2"T1 
|3"T1 13"2l [53"] . Then, protons or deuterons in a definite eigen state should 
remain in the same state in the course of time, while heavy atoms oscillate 
slowly, in an adiabatic hyperpotential depending on the proton/deuteron 
state, through the coupling term. This separation is relevant for KHCO3 
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and KDCO3 because adiabatic potentials for different proton or deuteron 
states do not intersect each other. In fact, the separation is rigorously 
exact in the ground state, since protons or deuterons should remain in 
this state for ever, if there is no external perturbation. Then, protons 
are bare fermions, while deuterons are bosons, so quantum correlations 
should be different for the two crystals. 



2.2 Macroscopic states 

Consider a crystal composed of very large numbers N a , N^, N c (TV = 
N a N^N c ) of unit cells labelled j,k,l, along crystal axes (a), (b), (c), re- 
spectively. The two dimers per unit cell are indexed as j, k, I and j', k, I, 
respectively, with j — j' . For centrosymmctric dimers, there is no per- 
manent dipolar interaction, so interdimer coupling terms and phonon dis- 
persion are negligible [22J, GUI [35] . The eigen states of the sublattice of 
protons can be therefore represented in a rather simple way with the basis 
sets of eigen states for isolated dimers 

A dimer (H1,H2) or (D1,D2) is modelled with coupled centrosym- 
metric collinear oscillators in three dimensions, along coordinates {aijki} 
and {a2jki} (a = x,y,z). The center of symmetry is at {aojki}- The 
mass-conserving symmetry coordinates independent of j, k,l, and their 
conjugated momenta, 



a s = —j= (ati - a 2 + 2a ), P sa = —= (P la - P 2a ), 
oc a = -^=(a 1 + a 2 ), P &a = -j= (Pi a + P 2a ), 



(2) 



lead to uncoupled oscillators at frequencies huj sa and ftw aa , respectively, 
each with m = 1 amu. The difference (huj sa — fno aa ) depends on the 
coupling term (say A Q ). The wave functions, {*njki( a a)}, {*n'jkl( a s - 
V2ao)}, cannot be factored into wave functions for individual particles, 
so there is no local information available for these entangled oscillators. 
Then, the wave functions of protons (fermions) or deuterons (bosons) are 
subject to the symmetrization postulate of quantum mechanics |36j . 



2.2.1 Protons 

The degenerate ground state of indistinguishable fermions must be anti- 
symmetrized with respect to site permutation. For this purpose, the wave 
function is written as a linear combination 
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Ojkli 



V2- 



II *0jklK) k jk iK - V2a ) ± *^ jkl (a s + V2a )l (3) 
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and the antisymmetrized state vectors with singlet-like \S) or triplet-like 
\T) spin symmetry are: 

\0jkl+)®\S) = |e 0j ki+> ® [| UU) - 1 4iT 2 >] ; 
\ojki-)®\T) = |e 0jkl _)®-i={|TiT2> + |4il2) (4) 



)=[\ Tiia> + UiT 2 )]}. 



V2 

The oscillators are now entangled in position, momentum, and spin. In 
contrast to magnetic systems [37], there is no level splitting, so the sym- 
metry related entanglement is energy-free. It is also independent of A a . 
Furthermore, as there is no significant exchange integral for protons sep- 
arated by ~ 2.2 A L 35], protons are not itinerant particles and there is no 
sizeable energy band width. 

Consider now the whole sublattice of protons. The spatial period- 
icity leads to collective dynamics and nonlocal observables in three di- 
mensions. With the vibrational wave function for the unit cell j, k, I, 
namely Sojkir = ©ojkir ± ©oj'kir, where r = "+" or "— " for singlet-like or 
triplet-like symmetry, respectively, trial phonon waves are written as 

1 n c N b N a 

Sor(k) = ^=^^^S 0jk lrCXp(ik-L), (5) 

1=1 k=l j=l 

where k is the wave vector, L = ja + kb + Ic is the lattice vector, a, 
b, c, are the unit cell vectors. This equation could represent collective 
dynamics of (H,H) dimers if they were composed bosons in a crystal made 
of indistinguishable dimer entities (KHC03)2. However, x-ray and neu- 
tron diffraction show that the structure is composed of monomer entities 
(KHCO3) related to each other through symmetry operations. The prob- 
ability density of each atom is equally distributed over all equivalent sites 
and, conversely, the probability density at each site includes contributions 
from all indistinguishable nuclei of the same kind. Consequently, antisym- 
metrization applied to the sublattice of indistinguishable fermions leads 
to 

k • L = modulo 2tt. (6) 

This means that there is no phonon (no elastic distortion) in the ground 
state. This symmetry-related "super-rigidity" [17] is totaly independent 
of any proton-proton interaction. The state vectors in three dimensions 
can be then written as: 

|H + ) = |So+(k = 0))®|5); () 
|H_> = |So_(k = 0)>®|T>. [n 

The wave functions So r (k = 0) represent collective oscillations of the 
super-rigid lattice as a whole, with respect to the center of mass of the 
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crystal. Finally, the ground state of the sublattice is a superposition state 

as: 

VAT|S 0+ (k = 0))®|S); 

vW|H _(k = 0))®|T). {0> 

This ground state is intrinsically steady against decoherence. Irradiation 
by plane waves (photons or neutrons) may single out some excited states. 
Entanglement in position and momentum is preserved, while the spin- 
symmetry and super-rigidity are destroyed. However, the spin-symmetry 
reappears automatically after decay to the ground state, presumably on 
the time-scale of proton dynamics. Consequently, disentanglement reaches 
a steady regime such that the amount of transitory disentangled states is 
determined by the ratio of density-of-states for the surroundings (atmo- 
sphere, external radiations...) and for the crystal, respectively. This ratio 
is so small that disentangled states are too few to be observed. Never- 
theless, they allow the super-rigid sublattice to be at thermal equilibrium 
with the surroundings, despite the lack of internal dynamics. 

The main source of disentanglement is the thermal population of ex- 
cited proton states. However, even at room temperature, the thermal 
population of the first excited state (< 1% for 7OH w 1000 cm -1 ) is of 
little impact to measurements. 

2.2.2 Deuterons 

For the sublattice of bosons in the isomorphic crystal of KDCO3, Eqs ([3]), 
Q and (|6|) are not relevant. There is neither spin-symmetry nor super- 
rigidity. Dynamics are represented with symmetry coordinates ([2]) leading 
to wave functions 

©ojki - J] *0jklK)* 8 0jkl (a 8 - V2a ), (9) 

a 

and phonon states 

, N c N b N a 

l Dk ±) = 7X7 E E E l S °iki±) exp(zk • L), (10) 

i = i k= i j=x 

where Sojki± = ©ojki ± Goyid- Needless to say, the H and D nuclei have 
the same number of degrees of freedom (12 per unit cell), but the sym- 
metrization postulate shrinks the size of the allowed Hilbert space from 
~ 12^ for bosons to ~ 12JV for fermions. 

3 Probing quantum entanglement with neu- 
trons 

Neutrons (spin 1/2) are unique to probing the spin-symmetry of macro- 
scopic proton states (JT]). However, quantum entanglement is extremely 
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fragile, as it is not stabilized by any energy. Only "noninvasive" experi- 
ments, free of measurement-induced decoherence, are appropriate [4]. For 
neutron scattering, this means (i) no energy transfer (ii) no spin-flip and 
(iii) particular values of the momentum transfer vector Q. (Q = k; — kf , 
where k; and kf are the initial and final wave vectors, respectively.) 




Figure 2: Diffraction patterns of KHC0 3 at 30 K (left) and 300 K (right) 
in the (a*,c*) reciprocal plane at k = 0, after Ref. [17] . The arrows 
emphasize ridges of intensity parallel to Q z and perpendicular to the dimer 
plane (dash lines along Q x )- 

Noninvasive diffraction events preserving the super-rigidity occur when 
components of the momentum transfer along proton coordinates (Q x > Qy, 
Q z ) match a node of the reciprocal sublattice of protons, so the momen- 
tum transfer is the same at all sites. Then, neutrons probe super-rigid 
states without any induced distortion. The only information conveyed by 
such events is the perfect periodicity of the sublattice, so the Debye- Waller 
factor is equal to unity at any temperature. In addition, thanks to the 
spin-symmetry, the scattered intensity is proportional to the total cross- 
section an ~ 82.0 b .17, 39 . Otherwise, if the matching condition is not 
realized, (i) the spin-symmetry is destroyed, (ii) the intensity scattered 
by protons is proportional to the coherent cross-section chc ~ 1-76 b, 
(iii) Bragg-peaks are depressed by the usual Debye- Waller factor for non- 
rigid lattices. Quantum correlations can be therefore clearly observed, 
thanks to the dramatic enhancement factor ch/chc * 45. Furthermore, 
the intensity scattered by heavy atoms, proportional to cr C Kco 3 ~ 27.7 b, 
is depressed by the Debye- Waller factor: exp — 2T / Kkco 3 (Q)- Therefore, 
the contribution at large Q-values is much weaker than that of the en- 
tangled sublattice and the contrast of intensity is enhanced at elevated 
temperatures. 
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The dashed lines in Fig. [T]show that proton sites are aligned along x 
and y, but not along z. Consequently, the noninvasive condition can be 
realized for Q x and Q y , whereas Q z does not coincide, strictly speaking, 
with a nod of the reciprocal lattice of protons. The diffraction pattern 
is composed of rods of diffuse scattering parallel to Q z , cigar-like shaped 
by the Debye- Waller factor. A detailed analysis of the scattering function 
reveals a complex pattern of rods of intensity, in full agreement with 
observations [T71 [T5] . For example, Fig. [5] shows slices of the diffraction 
patterns of KHCO3 at 30 and 300 K, in the (a*,c*) reciprocal plane, at 
Q y = (k = Qy/b* = 0). The cigar- like shaped rods are observed at 
Q x = and ±(10.00±0.25) A -1 , in accordance with the distance of « 0.6 
A between double lines of proton. These rods are evidences of macroscopic 
quantum correlations in the sublattice of protons. As anticipated from 
(fTQ|) . these rods are not observed for KDCO3 [16] . 




Figure 3: Projections on the (a, b) plane of the L (bottom) and R (top) 
configurations. 

At low temperatures, the rods are partially hidden by the anisotropic 
continuum, centered at Q = 0, due to incoherent scattering. At 300 
K, this continuum is dramatically depressed by the Debye- Waller factor, 
while diffraction by the super-rigid sublattice is largely unaffected. Quite 
remarkably, quantum correlations are not destroyed by interconversion. 
The sublattice is therefore a superposition of macroscopic states, \L) and 
\R), corresponding to the structures sketched in Fig. [3] 

Theory and experiments show that protons (deuterons) in the crystal 



10 



field are not individual particles possessing properties in their own rights. 
Instead of that, they merge into macroscopic states. The theory suggests 
that such matter fields, intrinsic to translational invariance, should occur 
in many crystals. However, the diffraction patterns in Fig. [2] are easy 
to interpret because protons are indistinguishable and hydrogen bonds 
are aligned (Fig. [1]). For more complex structures, seeking rods of dif- 
fuse scattering in diffraction patterns could be problematic. In the next 
section, we propose a complementary approach to evidence the nonlocal 
nature of protons and deuterons. 

4 Isotope mixtures 

In the present section we examine whether macroscopic states exist in 
isotope mixtures KH p D 1 _pC0 3 . We anticipate that any superposition of 
d?]) and (fTU)) should preserve the P2\/a space group symmetry of the crys- 
tal, whereas any statistical distribution of isotopes, here or there, should 
destroy the translational invariance. Because neutrons are scattered by 
nuclei, the diffraction pattern should discriminate between these two pos- 
sibilities. By contrast, the P2±/a symmetry evidenced by x-ray would not 
exclude, per se, a random distribution of isotopes, for they have identical 
electronic structures. 

A nearly cubic specimen (3x3x3 mm 3 ) was cut from a large crystal 
grown by cooling slowly a saturated solution in a mixture of H2O (sa 75%) 
and D 2 (« 25%). Measurements were conducted with the four-circle 
diffractometer 5C2 at the Orphee reactor (Laboratoire Leon-Brillouin) 
[40j . The crystal wrapped in aluminum was maintained at (300 ± 1) K. 
Data reduction was carried out with CRYSTALS [HI 132"] . 

Table 1: Unit cell parameters in A and units determined with single- 
crystal neutron diffraction at 300 K. A = 0.8305 A, space group P2i/a. 
The variance for the last digit is given in parentheses. 



A preliminary inspection of intensities at the absent reflections (OfcO for 
k 7^ 2n and hOl for h ^ 2n) confirms the P2\/a space group assignment. 
The unit cell parameters in Table Q] are identical to those of KHCO3. 
A complete determination of the positional and thermal parameters is 
currently in progress. 



KH . 75 D . 25 CO 3 KHCO3 [17 



6(A) 
c(A) 




15.180(1) 15.180(1) 
5.620(1) 5.620(4) 
3.710(1) 3.710(4) 
104.67(1)° 104.67(5)° 
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The P2i/a symmetry is at variance with a distribution of H and D 
nuclei at random. It is consistent with the existence of four sublattices, 
Hl, H#, Dl, Dfl, analogous to those sketched in Fig. [31 such as: 

RD(p,p H ,p D ,)= . , 

p[(l-p H )E L + PH E R \ + {l-p)[(l-p D )D L +pD'D R \. y > 

With a diffractometer based at a reactor source one measures a lim- 
ited range of the reciprocal space, typically |Q| ^ 5 A -1 . The matching 
condition for noninvasive measurements is realized only at Q x = Q y = 0, 
so for almost all measured Bragg-peaks, the spin symmetry is destroyed, 
the cross-section is er Hc , an d the sublattice of protons behaves exactly 
the same as a non-rigid lattice of bosons (fTU)) . The four sublattices are 
therefore consistent with a superposition of separable macroscopic states 
such as: |Hk±)z,, |Hk±)_R, |Dk±)i, |Dk±)fl- There is no direct information 
as to whether quantum correlations for proton states hold in the mixed 
sublattice. This point deserves further investigations. 

5 Conclusion 

The P2x/a space group symmetry of the KH p Di_ p C03 (p « 0.75) crys- 
tal determined with neutron diffraction shows that H and D nuclei are 
not individual particles located at definite sites and mutually exclusive. 
Wave-like protons and deuterons are equally distributed over all indis- 
tinguishable sites and, conversely, each site is occupied by all wave-like 
H and D nuclei with total occupancies p and (1 — p), respectively. This 
nonlocal behaviour allows the sublattice H p Di_ p to match the symme- 
try of the electronic structure that is unaffected by isotope substitution. 
The superposition of macroscopic states sounds as weird as a Cat in a 
superposition of |alive) and |dead) states. 

To the question "Where are protons and deuterons?", our answer is 
that the particle-like representation, implicitly meant by "Where" , is ir- 
relevant. Only the wave-like, or matter-field, representation is permitted 
by the crystal structure. The theory is based on fundamental principles of 
quantum mechanics and it is free of any ad hoc hypothesis or parameter. 
The nonlocal nature of protons and deuterons could be therefore of gen- 
eral significance. Neutron diffraction by hydrogen bonded crystals with 
partial isotope substitution could be a well suited and powerful approach 
to study nonlocal behaviors. 
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